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^ ■ Abstract 

P ■ 

The dynamics of two planar elastic pendula mounted on the horizontally excited platform have been studied. We give 
evidence that the pendula can exhibit synchronous oscillatory and rotation motion and show that stable in-phase 
and anti-phase synchronous states always co-exist. The complete bifurcational scenario leading from synchronous 
to asynchronous motion is shown. We argue that our results are robust as they exist in the wide range of the system 
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parameters. 
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1 Introduction 

The elastic pendulum is a simple mechanical system which comprises heavy mass suspended from a fixed point by a 
light spring which can stretch but not bend when moving in the gravitational field. The state of the system is given 
by three (spherical elastic pendulum) or two (planar elastic pendulum) coordinates of the mass, i.e. the system has 
■ three (spherical case) or two (planar case) degrees of freedom. The equations of motion are easy to write but, in 
general, impossible to solve analytically, even in the Hamiltonian case. The elastic pendulum exhibits a wide and 
surprising range of highly complex dynamic phenomena [U[2l[3l|4l[5l[6l[Zl[8l[9l [10l [Til [T2J [13l [14l [15l [16j . 

For small amplitudes perturbation techniques can be applied, the system is integrable and approximate analytical 
solutions can be found. The first known study of the elastic pendulum was made by Vitt and Gorelik [17] . They 
considered small oscillations of the planar pendulum and identified the linear normal modes of two distinct types, 
vertical or springing oscillations in which the elasticity is the restoring force and quasi-horizontal swinging oscillations 
in which the system acts like a pendulum. When the frequency of the springing and swinging modes are in the 
ratio 2 : 1, an interesting non- linear phenomenon occurs, in which the energy is transferred periodically back and 
forth between the springing and swinging motions [U [2J 02 SJ [6] . The most detailed treatment of small amplitude 
oscillations of both plane and spherical elastic pendula is presented in the works of Lynch and his collaborators 
[TUl HU . For large finite amplitudes the system exhibits different dynamical bifurcations and can show chaotic 
behavior E M M H3 M ■ 

The dynamics of elastic pendulum attached to linear forced oscillator has been studied by Sado [18] . She has 
shown a one parameter bifurcation diagrams showing different behaviour of the systems (periodic, quasiperiodic 
and chaotic). According to our knowledge this is the only study of considered systems, but one can find a lot 
of papers concerning dynamics of classical pendulum attached to linear or non-linear oscillator. Hatwal et al. 
[T51 |2T)I |2"T] gives approximate solutions in the primary parametric instability zone, which allows calculation of the 
separate regions of periodic solutions. Further analysis enables us to understand the dynamics around primary and 
secondary resonances [22, 23, 24, 25, 26 . Then the analysis was extended to systems with non-linear base where non- 
linearity is usually introduced by changing the linear spring into nonlinear one [27l [28l [29J [3Q| or magnetorheological 
damper [31] . Recently the complete bifurcation diagram of oscillating and rotating solutions has been presented 
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Figure 1: Model of the system 



[32] • Dynamics of two coupled single- well Duffing oscillators forced by the common signal has been investigated 
in our previous papers |33l 134) . We have shown the detailed analysis of synchronization phenomena and compare 
different methods of synchronization detection. 

In this paper we study the dynamics of two planar elastic pendula mounted on the horizontally excited platform. 
Our aim is to identify the possible synchronous states of two pendula. We give evidence that the pendula can 
synchronize both in the oscillatory and rotational motion moreover in-phase and anti-phase synchronizations co- 
exist. Our calculations have been performed using software Auto-07p [35J developed for numerical continuation of 
the periodic solutions and verified by the direct integration of the equations of motion. We argue that our results 
are robust as they exist in the wide range of the system parameters. 

The paper is organized as follows. Sec. 2 describes the considered model. We derive the equations of motion 
and identify the possible synchronization states. In Sec. 3 we study the stability of different types of synchronous 
motion. Finally Sec. 4 summarizes our results. 

2 Model of the system 

The analyzed system is shown in Fig. [TJ It consists of two identical elastic pendula of length l , spring stiffness ki 
and masses m, which are suspended on the oscillator. The oscillator consists of a bar, suspended on linear spring 
with stiffness k\ and linear viscous damper with damping coefficient c\. The system has five degrees of freedom. 
Mass M is constrained to move only in vertical direction and thus is described by the coordinate y. The motion 
of the first pendulum is described by angular displacement ip and its mass by coordinate X2, that represent the 
elongation of the elastic pendulum. Similarly the second pendulum is described by angular displacement (f> and its 
mass by coordinate X3. Both pendula are damped by torques with identical damping coefficient C2, that depend 
on their angular velocities (not shown in Fig. [!}. The small damping, with damping coefficient C3 is also taken for 
pendula masses. The system is forced parametrically by vertically applied force F(t) = Fq cos j4, acting on the bar 
of mass M, that connects the pendula. Force Fq denotes the amplitude of excitation and v the excitation frequency. 
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The equations of motion can be derived using Lagrange equations of the second type. The kinetic energy T, 
potential energy V and Rayleigh dissipation D are given respectively by: 

T = \{M + 2m)y 2 + \mx\ + \m{h + y ws ti + x 3 ) 2 <ft 2 + myx 3 cos0 - my<fi(l + y ws ti + x 3 ) sin</>+ (1) 
+ ^mxl + ^m(l + y ws t2 + x 2 ) 2 <j> 2 + myx 2 costp - my<p(l + Vwsti + x 2 ) sintp 

V = -mg(l + y ws t2 + x 2 )cosf- mg(l + y ws t2 + x 3 ) cos </> + mg(l + y wst2 ) + mg(l + y ws t2)+ (2) 
+ ^ fc i(y + ywstif + ^k 2 (y wst2 + x 2 ) 2 + ^k 2 (y wst2 + x 3 ) 2 - (M + 2m)gy 

D = l -C 2 ^ 2 + \c 2 i? + \c 3 x\ + \c 3 x 2 (3) 

where c 3 is the damping coefficient of the pendulum mass and y wst i — ( M +^ m )g ; y wst2 = ^ represent static 
deflation of mass M and pendulums' mass m respectively. The system is described by five second order differential 
equations given in the following form: 

m(l + Vwst2 + x 2 ) 2 ip + 2m(l + y wst2 + x 2 )(px 2 - my(l + y wst 2 + x 2 ) sin ip+ (4) 
+mg(l a + y wst2 + x 2 ) simp + C 2 ip = 

m(lo + Vwst2 + x 3 ) 2 4> + 2m(l + y wst2 + x 3 )4>x 3 - my(l + y wst2 + x 3 ) sin 0+ (5) 
+mg(l + y wst 2 + x 3 ) sin</> + C 2 4> = 

mx 3 + mycos<p- m<fi 2 (l + y ws t2 + x 3 ) - mg cos <j) + k 2 (y wst2 + x 3 ) + C 3 x 3 = (6) 

mx 2 + mycosip - mip 2 (l + y wst2 + x 2 ) - mgcosp + k 2 (y wst2 + x 2 ) + C 3 x 2 = (7) 

(M + 2m)y + mx 3 cos0 - 2mx 3 4>sm(j) - m(l + y ws t2 + x 3 )(j)sm(j)- m(l + y ws t2 + x 3 )4> 2 cos0+ (8) 

+mx 2 cos Lp 2mx 2 ipsmip - m(l + y wst 2 + x 2 )ip sinp - m(l + y ws t2 + x 2 )ip 2 cosip+ 

— (M + 2m)g + ki(y + y wst i) + C x y - F cosz/i = 

In the numerical calculations we use the following values of parameters: M = 10 [kg], m — 0.2 [kg], lo — 0.24849 [m], 
fci - 1642.0 [£], k 2 = 19.7 [^], ci = 13.1 [^] , c 2 = 0.00776 [Nms], c 3 - 0.49 [^], y wM = 0.062 [m], y wst2 - 
0.1 [m]. 

Introducing dimensionless time r = w\t, where lo 2 = M ^ 2m 1S tne na tural frequency of mass M with the attached 
pendula, we obtain dimensionless equations of motion written as: 

* + 77-^ ^X2-j-^ T 7sin* + - SlnVl/ ^ ^* = (9) 

(l+y2st+X2) (l+y2st+X2) (1 + V2st + XV (l+y2st+X2Y 
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i 2/3 2 i . 8? .. . x sin$ a 2 i 

— 7$>X3-t, — 77sin$ + - r + - ^$ = 10 

(i + y2 S t + X3j (i + y2 S t + X3j (i + y2 S t + X3) (l + yist + xz) 

X3 + StCOS $ - i±i^±li$2 _ 1 cos $ + yst2 +X3+ a ^ = (n) 

P2 P2 P2 

X2 + %i cos * - 1 + ?fet+X2 vj/2 _ 1 cog ^ + ^ + ^ + ^ = Q (12) 

P2 P2 "2 

1 + 4?X3 cos $ - %^X3* sin $ - ( 1 +W'«+ X3)a $ sin $ - a+w.t+X3)<> $2 cos $ + fc 2 cos ^ + (13 ) 

Pi Pi Pi Pi Pi 

W~ X2*sin* -J *sm* * cos* - -5- +7 + yut + «i7 - q cos fir = 

Pi Pi Pi Pi 

where uj 2 - — , w 4 - ^, fj, - — , Pi - — , p 2 - — , a - jjip^, q - u i ll(M+2m) , "l - wl(M +2m)> a2 ~ ^ij' 

03 m/owl ' yist lo ' ^ 2st lo ' 7 Jo' 7 '0^4' ' 7o^J' X3 Jo' X3 Jo"2' * 3 Zo^J'* 2 Jo' 

X2 = T^ 2 -, X2 = 7%, * = V, * = A * = 4, $ = 0, $ = ^, 3> = 4 

The dimensionless parameters of the system have the following values: 8\ = 2, 82 = 1-58, ct\ = 0.1, 0:2 = 0.01, 
a 3 = 0.1, a = 0.0192, y lst = 0.25, y 2s t = 0.4. 

We study system (9-13) in order to detect possible synchronization ranges. There are two basic types of 
synchronous motion, which are depicted in Fig. 2(a,b). The pendula can synchronize either in-phase or in anti- 
phase with each other, i.e., f = 9 or f = — 0. In both mentioned cases the forces acting in vertical direction 
on mass M are identical (there are no forces in horizontal direction), hence the energy transmitted between mass 
M and pendula in in-phase and anti-phase motion is also identical. If there is an in-phase synchronization, the 
anti-phase also coexists in the same range of parameters. The accessibility of in-phase and anti-phase motion is 
governed only by initial conditions. The pendula's masses are always synchronized in the in-phase with each other, 
i.e., \2 = X3- The anti-phase configuration of the masses is not observed (X2 = — X3) with the oscillating pendula. 
The anti-phase synchronization of masses is possible when the pendula are in equilibrium positions, then the sum 
of forces transmitted to mass M is equal to zero. 
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Figure 2: Possible synchronization (a) in-phase, (b) in anti-phase 

3 Stability of synchronous motion 

3.1 Synchronous solutions in two dimensional parameters space 

In this section we study the stability of the observed synchronous oscillations and rotations of the pendula. We 
present the bifurcation diagrams calculated in two-parameter space: amplitude q versus frequency /i of excitation. 
We focus our attention on determining the regions of synchronous stable motion and bifurcations that lead to its 
destabilization. We consider the state of the system in the following range q £ [0.0, 1.2] of forcing amplitudes 
and frequency of excitation belonging to the range /i £ [0.3, 1.2], which cover the possible resonances in the 
system. Resonance should be observed when the frequency of excitation comes close to the natural frequencies: 
of mass M equal to /x^f = 1, of pendula fi p — 0.50 and pendulum mass /x pm = 0.79. We describe synchronous 
solutions with respect to the forcing period according to ratio r : 1 — s : 1, where r and s denotes number of 
forcing periods, for which pendula and pendula masses perform period one motion. Fig. [3] presents two parameter 
bifurcation diagram, obtained by direct integration of (9-13). It shows the existence of synchronous, asynchronous 
motion and equilibrium solutions. As soon as we have a lot of coexisting solutions to hold clearance of Fig. [3] we 
do not distinguish which type of synchronous or asynchronous we find. By synchronous solution we mean, that 
both pendula are in complete synchronization state, i.e., their amplitudes and frequencies are identical. For low 
amplitudes of excitation, the only solution is equilibrium, which turns into synchronous or asynchronous solution 
as the frequency of excitation increases. The detailed analysis of synchronous solutions (shaded in grey color) was 
performed using continuation software Auto-07p [35] ■ We calculate the stability borders of each identified case, i.e., 
the ranges inside which the given motion is stable. The first periodic solution is observed for frequency of excitation 
equal to /i = 0.406 and for amplitudes of excitation above q = 0.709. This periodic solution is shown in Fig. |4ja) is 
identified as synchronous oscillations of pendula and pendula masses locked 1:1 — 1:1 with forcing. This solution 
is destabilized by saddle-node (green line), period doubling (blue line) and Neimark-S acker (red line) bifurcations 
curves. The continuation reveals that for small range of parameters, around the frequency of excitation close to 
the natural frequency of pendula, this solution coexists with synchronous 2:1 — 1:1 oscillations. Synchronous 
oscillations 2:1 — 1:1 are destabilized by saddle-node bifurcation curve then by Neimark-Sacker and pitchfork 
symmetry breaking (SB2) bifurcations. In the investigated system we distinguish two different symmetry breaking 
pitchfork bifurcations one of them (SB2) brokes symmetry between the pendula, the second one (SB1) brokes 
the symmetry of each pendula but their motion remains identical (36]. As the frequency of excitation increases 
we observe either asynchronous motion or equilibrium. With further increase of excitation frequency we observe 
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asynchronous behavior, which change into two small regions of synchronous rotations 3:1 — 1:1. We show it in 
Fig. HJf) and this area is bounded by Neimark-Sacker, period doubling and saddle-node bifurcations. This solution 
coexists with synchronous 2:1 — 1:1 rotations, presented also in Fig. @Jf). The stability region for this solution 
is bounded by pitchfork SB1 bifurcation from the left and right, Neimark-Sacker from above, and saddle-node and 
Neimark-Sacker bifurcations from the right. Both these solutions coexist in small range of considered parameters 
with another synchronous rotations 4 : 1 — 1 : 1, presented in Fig. 0Jf). The synchronous motion destabilizes 
from the right by saddle-node and pitchfork SB2 curves, from above by Neimark-Sacker, and from the left by 
Neimark-Sacker, saddle-node and pitchfork SB2 curves. 

Around /i w 0.8, where the resonance of pendulum masses occur, the system possesses rich dynamics, which 
results in the coexistence of different synchronous together with asynchronous solutions. This includes synchronous 
rotations 2:1 — 1:1 depicted in Fig. |U[c,d) and synchronous 3:1 — 3:1 rotations of pendula and pendula masses 
presented in Fig. B|e). The third which was found solution is synchronous half-rotations 1:1 — 1:1 (Fig. |4jh)), 
for which both pendula stop before approaching stable and unstable equilibrium transferring the whole energy into 
displacement of pendula masses. This multistability causes that it is hard to compare the bifurcation diagrams 
from the direct integration and Auto-07p. In the case of rotations 2:1 — 1:1 the synchronous motion is destabilzed 
from the right by saddle-node and pitchfork SB2 curves, from above by Neimark-Sacker, and from the left by 
saddle-node, Neimark-Sacker and pitchfork SB2 bifurcations. Synchronous rotations 1:1 — 1:1 loose stability by 
pitchfork SB2 from the right, and by Neimark-Sacker and period doubling from the left. The synchronous rotations 
3:1 — 3:1 are mainly destabilized by pitchfork SB2 from the right and by pitchfork SB2 and period doubling 
from the bottom, and by period doubling, pitchfork SB2, saddle-node and Neimark-Sacker from the left. From this 
solution, through period-doubling bifurcation we find synchronized rotations 6:1 — 6:1, shown in Fig. [4je). This 
solution is destabilized from above by period-doubling bifurcation, from the left through pitchfork SB2 bifurcation, 
and from below through saddle-node bifurcation (not visible, since coincides with period-doubling boundary for 
rotations 3:1 — 3:1). 

As we pass through the resonance frequency of mass M equal to fi — 1, for higher amplitudes of excitation 
the only synchronous solution includes 2:1 — 1:1 synchronous rotations depicted in Fig. IHc,d). After the 
resonance, for amplitudes of excitation above q = 0.141, only asynchronous solutions are observed. Below this 
value, many small synchronous regions were found. This includes synchronous 1:1 — 1:1 oscillations and two 
regions of synchronous 2:1 — 1:1 oscillations, together with two regions of synchronous 2:1 — 1:1 rotations. The 
region of 1 : 1 — 1 : 1 oscillations is enclosed by saddle-node, Neimark-Sacker and pitchfork SB2 bifurcation curves. 
Oscillatory 2:1 — 1:1 motion destabilizes through pitchfork SB1 from above and Neimark-Sacker curves from 
below. This solution coexists for small range of parameters with 2:1 — 1:1 rotations, which motion is destabilized 
by period doubling and Neimark-Sacker from the left, and from the right by pitchfork SB2, Neimark-Sacker and 
period doubling curves. We observe the excellent correlation in these regions between the results from numerical 
continuation and direct integration. 
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Figure 3: (color online) The synchronous (black dots), asynchronous (red dots) and equilibrium (small grey crosses) 
solutions of system (9-13) in two parameters space: /i frequency and q amplitude of excitation. We calculate this 
plot by direct integration using 4th order Runge-Kutta algorithm. In rectangles (a-f) we highlighted regions of 
synchronous motion calculated in Auto-07p (see Fig. H]). 



3.2 One parameter continuation 

In this subsection we show one parameter continuation of four periodic solutions (two oscillating and two rotational) 
for fixed amplitude of excitation, as a bifurcation parameter we choose the frequency of excitation //. We start each 
path- following on the periodic solution and continue in two directions (forward and backward). In Fig. [SJJa-d) we 
present the synchronized oscillating periodic solutions, their regions of stability are shown in Fig. 0] (a). System @ 
II 3[) is given by five second order ODEs, hence the phase space is ten dimensional and at least five figures (amplitude 
of each degree of freedom) are necessary to show its complete dynamics. To decrease it we focus on the dynamics of 
first pendula (second pendula has the same amplitude in the synchronized state) and mass M. The first presented 
branch of periodic solutions in Fig. [5](a,b) is synchronous 2:1 — 1:1 oscillations, in previous subsection we 
show that this family is destabilized by Neimark-Sacker bifurcation from the right and from the left by pitchfork 
symmetry braking SB2. Switching the branch at pitchfork bifurcation enables us to find another stable branch of 
asynchronized oscillations 2:1 — 1:1, that looses its stability through the saddle-node bifurcation. After pitchfork 
symmetry breaking SB2 bifurcation the solutions of one pendulum is located at upper branch (see Fig. [SJb)) and 
the second pendulum on lower branch or vice- versa. The branch synchronous oscillations 2:1 — 1:1, shown in Fig. 
[5] (c,d), present much richer scenario than other. These oscillations destabilize from both sides through pitchfork 
SB2 bifurcation. When we switch branch in left SB2 point, we find family of stable asynchronized oscillations 
2:1 — 1:1. Finally, when the amplitudes of pendula reach zero their motion stops. When we continue in the 
opposite direction the stability is lost in pitchfork SB2 bifurcation. Another change of branch allows us to observe 
another asynchronous periodic solutions, for which first pendulum oscillates 2:1 — 1:1, second pendulum is 
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Figure 4: (color) Stable ranges of synchronous motion calculated in Auto-07p (see rectangles in Fig. [3]). Color 
of lines stand for different types of bifurcation: Neimark-Sacker (red), saddle-node (green), pitchfork SB1 (violet), 
pitchfork SB2 (yellow) and period doubling (blue). In region inside lines synchronous solutions are periodic and 
stable. 
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at rest (not shown here) and pendula masses oscillate 1 : 1 — 1 : 1 in asynchronized manner. One end of this 
stable branch destabilizes through saddle-node bifurcation and the second one by pitchfork SB2 bifurcation. As the 
frequency of excitation increases the stability of this solutions is regained through pitchfork SB2 bifurcation and 
lost again through saddle-node bifurcation. Note that for the mass M, the bifurcation points that are responsible 
for the destabilization of periodic solutions for synchronized oscillations 2:1 — 1:1 and asynchronized oscillations 
2:1 — 1:1, are placed very close to each other. When we switch the branch in the right SB2 bifurcation point of 
the synchronized oscillations 2:1 — 1:1, we find asynchronized solution of oscillations 2:1 — 1:1 that persists for 
small interval of excitation frequency. It destabilizes from above and below through period-doubling bifurcations. 
Switching the branch in lower period doubling bifurcation point enables us to observe asynchronized oscillations 
4 : 1 — 2 : 1, that destabilize through Neimark-Sacker bifurcation. The bifurcation diagram, shown in Fig. EJa-b), 
shows 3:1 — 3:1 rotational periodic solutions for q — 0.654. In contrary to the previous cases, to hold a physical 
meaning, on the horizontal axis we plot the amplitude of velocity. The stability region of this branch is shown 
in Fig. [He). This family of solutions looses its stability through period-doubling and pitchfork SB2 bifurcations. 
Switching the branch in the period-doubling bifurcation point, allows to observe another period doubling, which 
leads us to two branches of synchronized 6:1 — 6:1 rotational solution. This branch looses its stability once again 
through period-doubling bifurcation. After another switch of branch in period doubling bifurcation point, we reach 
two synchronized periodic 12:1 — 12:1 rotational solutions. They are stable in very narrow range of excitation 
frequency and loosing stability via Neimark-Sacker bifurcations. After switching the branch in right pitchfork SB2 
bifurcation point, we observe asynchronized rotations 3:1 — 3:1, that are stable in very small interval, finally 
loosing its stability through saddle-node bifurcation. In Fig. [6jc-d) we present synchronized rotations 4:1 — 1:1, 
that loose stability through pitchfork symmetry braking from the right and left. Switching the branch in both SB2 
bifurcation points let us to find asynchronous rotations 4:1 — 4:1, that are stable in very narrow interval of 
excitation frequency, loosing finally stability through saddle-node bifurcation. 

4 Conclusions 

In the system of two planar elastic pendula suspended on the excited linear oscillator one can observe both in-phase 
and anti-phase synchronization of the elastic pendula. In-phase and anti-phase synchronous states always co-exist. 
Pendula can synchronize during the oscillatory and rotational motion but only when their behaviour is periodic. 
We have not observed the synchronization of the chaotically behaving pendula. This result is contrary to the 
great number of chaos synchronization examples [37l EHJ [39] but confirms the results obtained in [40] where it has 
been shown that the forced Duffing's oscillators mounted to the elastic beam can synchronize only after motion 
become periodic. The synchronization of the chaotic motion of the pendula is impossible as the excited oscillator 
transfers the same signal to both pendula which cannot differently modify the pendula's motion. We also have 
not observed in-phase or anti-phase synchronization of the pendula when masses mi and 7713 are in anti-phase. 
With one parameter bifurcation diagrams we present a bifurcational scenario of synchronous solutions. We show 
a route from synchronous via asynchronous periodic solutions to quasiperiodic and chaotic behaviour. In this case 
the pendula in-phase or anti-phase synchronization is impossible as the pendula have different and non-constant 
lengths. In this case one can expect some kind of generalized synchronization but this problem will be addressed 
elsewhere |41j . 

We show two dimensional bifurcation diagrams with the most representative periodic solutions in the considered 
system. In the neighbourhood of the linear resonances of subsystems we have rich dynamics with both periodic 
and chaotic attractors[42j. Our results are robust as they exit in the wide range of system parameters, especially 
two dimensional bifurcation diagram can be used as a scheme of bifurcations in the class of systems similar to 
investigated in this paper. 
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Figure 5: 1 parameter continuation of fully synchronized: 2:1-1:1 oscillations ((a) mass M, (b) pendulum 1. q = 
0.899, /i = 0.455), 2:1-1:1 oscillations ((c) mass M, (d) pendulum 1, q = 0.654, n = 0.5). The continuous and 
dashed lines correspond to stable and unstable periodic solutions respectively. Abbreviations depicted following 
bifurcations: NS (Neimark-Sacker), PD (period doubling), SB1 (pitchfork SB1) and SB2 (pitchfork SB2). Other 
changes of the stability take place through the saddle-node bifurcations. 
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Figure 6: 1 parameter continuation of fully synchronized: 3:1 — 3:1 rotations ((a) mass M, (b) pendulum 1, 
q = 0.654, /i = 0.8), 4:1 — 1:1 rotations ((c) mass AI, (e) pendulum 1, q = 0.8, [i = 0.7). The continuous and 
dashed lines correspond to stable and unstable periodic solutions respectively. Abbreviations depicted following 
bifurcations: NS (Neimark-Sacker), PD (period doubling), SB1 (pitchfork SB1) and SB2 (pitchfork SB2). Other 
changes of the stability take place through the saddle-node bifurcations. 
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